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1 Introduction 



Massless charged particles do not seem to exist in nature, the Hghtest known charged particle 
being the electron. Apparently there is no fundamental theoretical principle that forbids their 
existence, although the consistency of relativistic quantum field theories describing massless 
charges and the related cancellation of infrared divergences are still open problems [1]. Surpris- 
ingly the consistency of those particles at the classical level is seemingly an open question and 
very little is known about their physical features: the fact that they travel with the velocity of 
light introduces mathematical, as well as physical, singularities that make the most standard 
basic instruments of Functional Analysis to solve Maxwell's equations - like the Green function 
method ~ difficult to apply. Basically what is known about classical massless charges is merely 
the electromagnetic field generated by a (free) charge in uniform linear motion, that is given by 
a (5-like shock-wave supported on a plane orthogonal to the trajectory, travelling at the speed of 



Here e is the charge of the particle, its four-velocity subject to the lightlike condition = 0, 
and (ux) stands for n^x^. 

In this paper we solve for the first time Maxwell's equations for a particle moving with the 
velocity of light along a generic accelerated trajectory and determine explicitly the electromag- 
netic field it produces. For the sake of simplicity we consider only bounded trajectories, but the 
new method we develop applies equally well to unbounded ones. Since the current density of a 
point-like particle is a Dirac (5-function, the solution makes sense only in the space of distribu- 
tions. The corresponding electromagnetic field J-'^'^, for which we give an analytical expression 
in (|4.10p . is Lorentz-covariant and respects causality and represents thus the generalization of 
the celebrated Lienard-Wiechert field from timelike to lightlike worldlines. The relevance of this 
result is first of all of conceptual character, since it proves that Maxwell's equations admit (exact) 
solutions in the space of distributions even if the worldlines are lightlike, although conventional 
methods fail in this case, as observed above. The knowledge of the exact field produced by 
an accelerated massless charged particle should eventually allow to settle the question whether 
or not a consistent classical radiation theory for such particles - preserving total energy and 
momentum - can be formulated; for preliminary analyses see e.g. The knowledge of 

this field represents furthermore the basis for a systematic analysis of possible interactions of 
this kind of particles, while the field (jl.ip describes only free particles. The results presented in 
this paper could thus shed new light on the open question of the possible existence of massless 
charged particles in nature, even though a definitive answer to this question can arise only in 
the framework of quantum field theory. Eventually the method developed in this paper to de- 
rive exact solutions for linear partial differential equations in the space of distributions can be 
applied also to other physically relevant problems, like the field produced by a massless dyon 
or the linearized gravitational field produced by a massless particle, like the photon. Detailed 
proofs and further developments will be presented elsewhere [9]. 

2 Heuristic analysis and the singularity surface F^^ 

A massless charged particle travels with the speed of light and, as long as its energy and mass 
are not directly involved in the physical process one considers, it appears natural to regard it as 
some limiting case of a particle that travels with speed V < 1. With this respect a speed V < 1 



light [His]: 




i^x" - u^xf 



(1.1) 
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in some sense plays a role complementary to a mass m > used frequently in quantum field 
theory to regularize infrared divergences. 

Timelike worldlines. Our starting point is therefore the standard Lienard-Wiechert field for 
a timelike worldline, that we write as 

pf.y ^ (J^.u ^ -^^,u ^ (2.1) 
where the Coulomb and radiation fields are given respectively by 



eU^ 

Att{ULY 
e 

4tt(ULY 



R'^ix) = -JZjhrVz mL)^ - {WLW] - (/i O z.). (2.3) 



Here we parametrized the particle's worldline Y^^ = Y^{X) with an arbitrary parameter A and 
defined L^'{x) = x-Y^',U^' = dYi'/dX, W^" = dU^'/dX. We adopt the convention [UL) = U^'L^ 
and similarly for iWL). In formulae ()2.2p and ()2.3p A is evaluated at its "retarded" value A(x), 
uniquely determined by the conditions 

l2 = (x-y(A))2 = 0, L° = - y°(A) > 0. (2.4) 

Notice that the square of the "four-velocity" is given by 

where V = \V\ = \dY /dt\ is the speed of the particle. The above fields satisfy the Maxwell 
equations 

in the sense of distributions : due to the presence of 5-like currents this is the only framework 
where these equations make sense. In particular the fields C^^ and R^^^ satisfy separately the 
Bianchi identities 

^[^.Rup]=^ = ^l^C,py (2.7) 

Lightlike worldlines. We consider now an arbitrary spatially bounded lightlike worldline 
= yi^i^X) and search for the causal solution of the associated Maxwell equations 

d[p^M = 0, d^T^'^ = ^ / ^^^"^ " ^^-^^ 

In analogy to the timelike worldline Y^ we define = dy^ /dX and = du^/dX where, since 
the wordline is lightlike, = 0. For future reference we introduce also a retarded parameter 
A*(x) defined analogously to <\2A\\ : 

/^(x) = x-y^(A*), /2 = 0, /°>0. (2.9) 



Loosely speaking we will derive the solution of equations (12. Sp considering a conveniently de- 
formed timelike worldline Y^{X) that in an appropriate limit - to be specified below - ap- 
proaches the given lightlike worldline y'^(A). In this limit we have X{x) — )• X*{x), U"^ ^ = 0, 
{UL) — )• [ul) and V{t) — )• v{t) = 1, where v{t) = dy/dt is the velocity of the lightlike worldline. 
For what concerns the fields (12. 2p and (12. Sp . since in this limit f/^ — t- the field C^^ appears 
to vanish, while R^'^ seems to admit a finite limit as it does not involve any power of U'^. Both 
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these conclusions are, however, of purely conjectural character due to the presence of zeros in the 
denominators of (|2.2p and (|2.c{p . These denominators are, in fact, powers of the scalar product 
(UL), that tends to (ul). But while (UL), given (j2.4p . vanishes only on the particle's worldline 
X = Y^{X), i.e. at fixed time t only at the particle's position, the scalar product {ul) vanishes 
along the physical singularity surface = Tf^{b, A) defined bjlj 

r'^(6,A) = y'^(A) + 6n'^(A), 6 > 0. (2.10) 

At fixed time t the spatial component of this surface is a Dirac-like string attached to the particle 
given by 

j{b,t) = y{t-b) + bv{t-b), b > 0. (2.11) 

The fact that (ul) vanishes on follows from the relations l'^ = {x — y(A))^ = 0, x'^ — y^{X) > 
and ti^(A) = 0. On the singularity surface F'^, that will play a crucial role in what follows, the 
denominators of (j2.2p and (j2.3p tend therefore to zero. This means that, if the above envisaged 
limits of C^^'^ and R'^'^ in some sense exist, these limits could involve a "special" tensor supported 
on r'^. Due to Lorentz-covariance, reparametrization invariance (A — t- A', 6 — )• {dX'/dX)b), 
antisymmetry in fi and ly and for dimensional reasons, this tensor can only be proportional to 

/OO roo 
dX biuf^w" -u''w^')6^ix-T{b,X))db. (2.12) 
OO Jo 

As a simple calculation shows this field entails furthermore another interesting property, that is 

d,P>'''=f, (2.13) 

where j'^ is the lightlike current defined in (j2.8p . For V ^ 1 the limit of the Coulomb field C^^'^ 
could thus only be proportional to P'^'^ since, as observed above, in the complement of it 
goes to zero. However, C'^'^ satisfies the Bianchi identity d^^C^p^ = and so also its (possible) 
limit kP^^'^ must satisfy it. But since, by inspection, dyP^p^ 7^ we conclude that k = 0: this 
implies that for a lightlike bounded worldline the Coulomb field must vanish identically. This 
argument does clearly not apply to the limit of the radiation field R'^'^, that could - and will 
indeed - develop a contribution proportional to P^^'^. 



3 A rigorous strategy 

To establish a well defined limiting procedure we have first of all to establish a feasible timelike 
regularized wordline Y^(X) for a given arbitrary lightlike bounded worldline y^{X). As will be 
clear in a moment, we need this regularized worldline to depend on a constant regularization 
parameter V. As simple as it may look, we define the regularized timelike wordline by 

Y\X)^^, Y{X)^y{X), 0<V<1. (3.1) 

The regularized orbit is thus the same as the orbit of y^(A), but the velocity of the particle is 
now V{t) = dY/dY^ = Vv{t) and hence its speed is a constant less than unity: \V{t)\ =V<1. 
The regularized wordline is thus timelike and the regularized field F'^'^ is given by (12. 2p and 
(j2.3p and depends on the constant parameter V. As this field is a distribution, it makes sense 
to analyse its limit as V tends to 1 in the sense of distributions. 

^For a uniform motion - which is necessarily unbounded - this surface collapses to the worldline of the particle 
and the following argument does not apply. In this case the Coulomb field admits actually the non-zero limit 
p.ip . For a qualitative description of the singularity curve (|2.11|) in a cyclotron motion see [S]. 
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Distributional limit. We want now to construct a solution of equations (j2.8|) in the distri- 
butional sense. In what follows we indicate the pointwise limit of a function with the standard 
symbol limy_^i, while we write the limit in the sense of (tempered) distributions as Limy_j.i. 
First of all we notice the trivial distributional limit, see equations (|2.6p and (|2.8p . 



Limv^i J^=f. (3.2) 

Since in the space of distributions derivatives are continuous operations, partial derivatives 
commute always with the distributional limit Limy^i. Therefore - and this is one of the key 
points of this paper - if the distributional limit of the field ()2.ip 



Limy_^i 

exists, applying the operation Limy_>.i to all equations in (12. 6p and interchanging the limits with 
the derivatives, one can conclude that the field J-^'^ satisfies automatically the Maxwell equations 
112. 8\) . In the rest of the paper we outline how the limit ()3.3p can be evaluated explicitly. 

We conclude this section noting that, thanks to (j2.5p . in the complement of F'^ the pointwise 
limit of (12.11) amounts to the function 



lim F^'^(x) = lim i?^^(x) = , l^'Uunw" - (wl)un - ^ u) = 7^^^(x). (3.4) 

4 The limiting procedure 

Limit of the Coulomb field. We discuss first the distributional limit of the Coulomb field 
Limy^iC'^'^ in ()2.2p . The evaluation of this limit requires to evaluate the ordinary limits 
limy_j.i C^*^ ((/?) = limy__i.i f C^^'^{x) ip{x) d'^x, where (p{x) is an arbitrary test function belong- 
ing to the Schwartz-space 5(M'^). For the reasons reported above these limits are conjectured to 
be all zero: using the techniques that we will develop below for the radiation field, an explicit 
calculation [9j shows indeed that these limits vanislo. We have thus 

Limy-.iC^'^ = 0. (4.1) 

Limit of the radiation field. In the remainder of this section we establish the existence of the 
distributional limit Limy_s.ii?'^'' and evaluate it explicitly, the result being (j4.10p . We perform 
the analysis explicitly for the electric component R^^ of the field; its magnetic component R^^ can 
be analyzed in an identical fashion. We must thus apply this field to a test function and consider 
the quantity R^^{(f) = J R^^{x) ip{x) d^x. To analyze this integral it is first of all convenient to 
disentangle the retarded parameter \{x) - transforming it in a "free" parameter - inserting the 
5-function identity 

2 j {UL)H{L^) 5{L^) d\ = j 5{\ - \{x)) dX = 1, L^' = x^" - y^(A), 

where H{ ■ ) stands for the Heaviside function. After that we perform the shift of variables 
j;'^ — )• + y^(A), integrate the (^-function over x^ and use then reparametrization invariance 
to choose as time variable t = Y^{X). Eventually we insert the regularized worldline (j3.ip and 
perform the rescaling t — t- t/V. In this way from (j2.3p we get 

RiO{^) = ^ J K\x)^{^+r,x + y)dtd^x, r = \x\, (4.2) 



^In this proof it is crucial that the trajectory is bounded. We stress that for a hnear uniform trajectory the 
limit (14.11) gives indeed the non vanishing result 
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= ^ : — —■ (4-3) 



where 

(a • x) (x* — rVv^) — r{r — V{v • x)) a* 
r (r — V{v ■ x)Y 

In expressions (j4.2p and (j4.3p the kinematical variables y, and a are respectively the position, 
velocity and acceleration of the lightlike worldline and are evaluated at the time t. We must now 
evaluate the limit of (j4.2p as V — )• 1. As long as y < 1 the denominator of (j4.3p vanishes only 
in r = 0, but for 1/ = 1 (at fixed t) it vanishes along the half-line x{b) = bv, b > 0, leading 
to a non integrable singularity in (j4.2p . This half-line is precisely the image of the singularity 
curve 7 in (j2.1ip . To isolate this line we perform in (j4.2p a final change of the space-coordinates 
X -H- (6, gi, (72); introducing at fixed t a basis Na orthogonal to -u = v{t): 

x = bv + q^N^, Na-N,3 = 6ais, Na-v = 0, a,/3=l,2. (4.4) 

In these coordinates the singularity curve 7 is represented simply by the conditions Qa = and 
6 > 0. Taking into account that v = 1, a ■ v = and that the Jacobian of the transformation 
()4.4p is 1, expression ()4.2p can be rewritten as 

i?*0((^) = ^ J {G' + W) if (^^ + r,bv + qaNa + y) dtdbcfq, (4.5) 



where r = ^^/W+q^, = Qi + Q2 and we introduced the functions 

= J^-V')ir^Vb)bY^ (4.6) 

(g2 + (1 - y2)52)2 

^ n^^jV^iVj gJ- g„(a.iV,)(fe-W)T;^ g'a^ . . 

r{r-Vb)^ r(r-F6)2 2r(g2 + (1 - ^2)^2) ' ^-^^ 



where 



J-J-a/3 = QaQp ^ <Jo/3- 



Under the limit V ^ 1 both integrals in (|4.5p converge now separately. The delicate point is 
represented by the (a priori non integrable) divergences arising at g = when V goes to 1: the 
interchange of the limit V ^ 1 with the integral over d^q must therefore be handled with care. 
In the G*-integral one can first perform the rescaling q^ — >■ — V^ qa and resort then to the 
dominated convergence theorem to move the limit V ^ 1 inside the g-integral. The resulting 
integral over d^q is then convergent, and elementary, and the test function if gets eventually 
evaluated at g = 0. For 6 < this integral is actually zero, because in this case under the above 
rescaling of qa the term (r + Vb) vanishes as y — )• 1. 

For what concerns instead the //^-integral, in the second and third terms of (j4.7p the limit 
can be taken trivially under the (7-integral sign. On the contrary the first term of (j4.7p is 
characterized by the singular behavior of its denominator, that for 6 > under V ^ 1 for 
small q behaves as l/(r — Vb)'^ — )• l/(r — 6)^ ~ l/(g2)2^ while its numerator is quadratic in q. 
For 1/ —7- 1 the corresponding integral seems thus logarithmically divergent. However, if before 
taking the limit V — )• 1 one performs first the integral over the polar angle of the g-plane 
and afterwards the integral over its modulus q = \/ qf + q2, the logarithmic divergence cancels 
because by symmetric integration the integral over ■(? of the matrix IIq,^ vanishes. In the so 
rearranged ff*-integral - enforcing again the dominated convergence theorem - one can take the 
limit y — )• 1 under the g-integral sign and the resulting integral is conditionally convergent in 



■^In the most simplest case an integral is said to be conditionally convergent, if it converges once one specifies 
a specific order of the integrations over the variables. 
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the sense just specified. In this way the hmit of (|4.5|) . that thanks to (|4.ip corresponds to the 
zO-component of the distributional limit l\'6.'6\i . becomes eventually the sum of two terms: 



OO JO 



r^{ip) = Urn R'\ip) = I dt I a'ip{t + b,bv + y)bdb (4.8) 




^ ] ^p(t + r,bv + qaNa+y) dtdbd^q, (4.9) 



r(r — 6)2 r(r — b) 2r 

where it is understood - we repeat - that the g-integral in (14.91) is conditionally convergent. 
The term in ()4.8p is supported on the singularity line (I2.11|) and by inspection one sees that it 
amounts precisely to |P*° ((/?), see (I2.12|) . To the distribution in (14. 9p one can formally apply 
all the coordinate transformations performed so far "backwards" , and in doing so one discovers 
that in the complement of it can be written as J TV'^{x) ip{x) d'^x, where the function TV'^{x) 
is defined in ()3.4p . This result is obviously not surprising, since in regions not containing 
singularities distributional and pointwise limits are equivalent. Repeating the same analysis for 
the magnetic components BJ'^ of ()2.3p . and taking into account that the above analysis holds for 
an arbitrary test function ip, one concludes that the limit (|3.3p exists and can be written as 



F^"" = i p/^^ + P(7l^^). (4.10) 

The "principal part" symbol V{TZ^'^) indicates the distribution defined in terms of the func- 
tion lZ^'^{x) given in (j3.4p . through the g- integration procedure around the singularity curve 7 
prescribed above. 



5 Analysis and comments 

We conclude the paper with an analysis of the most peculiar and salient features of the expression 
F^" in ([CT]1 . 

By construction F^^^ represents a distribution satisfying the Maxwell equations (12. Sp and 
is a Lorentz-covariant tensor. This last property is not at all trivial since our regularization 
()3.ip breaks Lorentz-invariance. There exist actually (under certain aspects) more complicated 
regularizations than ()3.ip . which have however the advantage of preserving manifest Lorentz- 
invariance [m [I2I [13] and which lead to the same result ()4.10p , see [9] . 

The Maxwell field J-^^^ is the sum of two non-regular distributions, i.e. distributions that are 
not represented by functions, the most striking feature being the appearance of the J-function 

The field J-^^ respects causality. For the contribution VilZ^^^) this follows from the fact that 
it is essentially the pointwise limit of a causal field, i.e. the Lienard-Wiechert field. The term 
1 pM!^ jg supported on the surface T^, that gives rise to the new physical singularity curve 7(6, t). 
From its explicit expression (j2.1ip one sees that the (measurable) velocity at time t of the point 
on the curve corresponding to the parameter b - given by the component of the vector d^{b, t)/dt 
orthogonal to the curve j{b,t) - is v{t — b): the singularity curve propagates therefore with the 
speed of light, in compatibility with causality. The causality character of the field (j4.10p is a 
priori non trivial, since it could not be derived relying on a causal Green function. 

The Hodge-dual of the distributional two-form P^'^ has a simple geometrical interpretation: 
it is the Poincare-dual of the surface r^(6. A), see e.g. [10]. The relation (I2.13P follows precisely 
from this property since the boundary of F^ is the worldline of the particle. 

Equation (I2.13P implies that in (I4.10p the electric fiux of the massless charge is carried half by 
the term ^ P^'^ and half by the principal part term, since from equations (|2.8p . (|2.13p and (|4.10p 
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it follows that d^jV{n^"') = \f. Notice, however, that while the field 2V{nf"') satisfies the 
second Maxwell equation in (|2.8|) . it would not satisfy the Bianchi identity, because d^^P^p^ ^ 
too. 

For the timelike Lienard-Wiechert field (j2.ip the electric flux - e - through a sphere centered 
at the (retarded) particle's position comes entirely from the Coulomb field, while the radiation 
field can be seen to give a vanishing contribution. On the contrary, once the distributional limit 
Limy^i on F'^'^ is taken, the Coulomb field disappears and the electric flux comes entirely from 
the "radiation fleld" ()4.10p . This represents a physically as well as mathematically interesting 
result: when taking the limit from a timelike to a lightlike wordline, the Gauss-law operation 
and the distributional limit are (maximally) non commuting operations. 

As last comment we observe that the energy-momentum tensor T^'^ associated to the Maxwell 
field (j4.10p . containing the square of a 5- function is ill-defined, as is the Poynting vector. This, 
however, does not imply that radiation theory for massless charges is necessarily inconsistent. 
In fact, even for timelike trajectories the tensor T^'^ - being the square of the standard Lienard- 
Wiechert field (12. ip - is not a distribution and has to be defined through a careful regulariza- 
tion/renormalization process, see [Tl]. It may be that such a process leads to a well-defined 
energy-momentum tensor also for lightlike trajectories. Given the explicit expression ()4.10p it 
should now be possible to provide for this question a definitive answer. 
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